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, We present two applications of quantum integrable systems. First, we predict that 

f***^ , it is possible to generate high harmonics from solid state devices by demonstrating 

^v^J i that the emission spectrum of a minimally coupled laser field of frequency u> to an 

impurity system of a quantum wire, contains multiples of the incoming frequency. 
Ch ' Second, by evaluating expressions for the conductance in the high temperature 

regime we show that the characteristic filling fractions of the Jain sequence, which 
occur in the fractional quantum Hall effect, can be obtained from quantum wires 
l^"} i which are described by minimal affine Toda field theories. 
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Dedicated to A. A. Belavin on his 60 th birthday 



1 Introduction 



In the context of conformal and massive integrable quantum field theories an im- 
pressive amount of non-perturbative techniques has been developed during the last 
25 years. Needless to say that the contributions of the Landau school has always 
been vital for that progress. The original motivation for considering such theories 
was to use the lower dimensional set up as a testing ground for general conceptual 
Oh. ideas and possibly to apply them in the context of string theory. As a consequence, 

most of the work in this area is often of a rather formal mathematical nature and 
lacks a link to direct physical application. So far this has not been a major issue, 
but lately the experimental techniques have advanced to such an extent that one 
can realistically hope to measure various physical quantities which can be predicted 
based on the developed approaches. 

Here we want to present two examples of such quantities. The first is concerned 
with the prediction of harmonic spectra when a three dimensional laser field is cou- 
pled to a one dimensional quantum wire™ To observe interesting phenomena in 
this context one needs to consider impurity systems and the concepts of integra- 
bility are so constraining that they lead one to consider mostly free systems with 
impurities. 

The second application we shall present is related to the computation of par- 
rel 

ticular values of the conductance of a quantum wire . There exist two established 
theoretical descriptions to compute the conductance, the Kubo formulcPEl which 
is the outcome of a dynamical linear-response theory and the Landauer-Biittinger 
theory^, which is a semi-classical transport theory. In both descriptions one can 
make full use of different non-perturbative techniques developed in the 1+1 dimen- 
sional quantum field theory context. In the first approach one of the key quantities 
involved is the curren t-curr ent correlation function which can be obtained from a 
form factor expansio n" ! ' ' ° > and in the second approach one requires the density 
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distributions which are accessible from a thermodynamic Bethe ansatz- (TBA) 
analysis. Here we will concentrate on the latter and show that the characteristic 
Jain filling fraction^, which occur in the fractional quantum Hall effeciP, can 
be obt aine d from quantum wires which are described by minimal affine Toda field 
theorieJEl 



2 Harmonic generation 

We commence by briefly explaining what har mon ics are. The first experimental 
evidence can be traced back to the early sixtiei^l Franken et al found that when 
hitting a crystalline quartz with a weak ultraviolet laser beam of frequency ui, it 
emits a frequency which is 2u>. Generalizing this phenomenon to higher multiples, 
one says nowadays that high harmonic generation is the non-linear response of a 
medium (a crystal, an atom, a gas, ...) to a laser field. The importance of harmonic 
generation is related to the fact that it allows to convert infrared input radiation 
of frequency to into light in the extreme ultraviolet regime whose frequencies are 
multiples of u> (even up to order ~ 1000, see e.g.LL^for a recent review). A typical 
experimental spectrum is presented in figure 1. 
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Figure 1. Harmonic spectrum for Neon for a Ti:Sa laser with A = 795nm. Measured at the Max 
Born Institut Berlin [15]. 
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In gases, composed of atoms or small molecules, this phenomenon is well-understood 
and, to some extent, even controllable in the sense that the frequency of the highest 
harmonic, the so-called "cut-off', visible in figure 1, can be tuned as well as the in- 
tensities of particular groups of harmonics. In more complex systems, however, for 
instance solids, or larger molecules, high-harmonic generation is still an open prob- 
lem. This is due to the fact that, until a few years ago, such systems were expected 
not to survive the strong laser fields one needs to produce such effects. However, 
nowadays, with the advent of ultrashort pulse s, th ere exist solid-state materials 
whose damage threshold is beyond the required^ intensities of 10 14 W/cm 2 . As 
a direct consequence, there is an increasing interest in such materials as potential 
sources for high-harmonics. In fact, several group s are currently inv esti gating this 
phenome non in systems such as thin crystals ^ 3^1 carbon nanotubes^l, or organic 
molecule ^ ! 21 ! 

We will therefore try to answer here the question, whether it is possible to 
generate harmonics from solid state devices and as a prototype of such a system 
we study a quantum wire coupled to a laser field. 

2. 1 Constraints from integrability 

In order to couple the laser field to the wire we need some seeds to attach the field 
and therefore we are naturally led to consider impurity systems. Essential quantities 
to compute are the transmission amplitudes. We commence by demonstrating how 
integrability puts severe constraints onto them. One exploits here one of the great 
advantages of integrability in 1+1 dimensional models, which is the well-known fact 
that the n-particle scattering matrix factorises into two-particle S-matrices, which 
can be determined by some constraining equations which are central to the entire 
subject, the Yang-BaxteJ 2 - ^ and boo tstrap equationfP^ Similar equations hold in 
the presence of a boundary^^SESI or a defecJ^EHI It is clear that with regard 
to the conductance a situation with a pure boundary, i.e. non-trivial effects on the 
constrictions, or purely transmitting defects will be rather uninteresting and we 
would like to consider the case when R and T are simultaneously non- vanishing. It 
will turn out that for that situation the Yang-Baxter equations are so constraining 
that not many integrable theories will be left to consider. Thus this section serves 
essentially to motivate the study of the free Fermion, which after all is very close 
to a realistic system of electrons propagating in quantum wires. 

We label now particle types by Latin and degrees of freedom of the impurity 
by Greek letters, the bulk scattering matrix by S, and the left/right reflection and 
transmission amplitudes of the defect by R/R and T/T, respectively. Then the 
transmission and reflection amplitudes are constrained by the "unitarity" relations 

R%m%(-0)+Ti£(6)f%(-9) = 5$5Z, (1) 

R£W%(-0) + Tg(0)B%(-9) = Q , (2) 

and the crossing-hermiticity relations 

R«(9) = R«(-6)* = S i3 (29)R?(iir - 9) , (3) 

T?(9)=f?(-9y=f?{in-6). (4) 
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The equations Q and (J2J also hold after performing a parity transformation, that 
is for R <-> R and T <-> f. 

Depending now on the choice of the initial asymptotic condition one can 
derive the following two non-equivalent sets of generalized Yang-Baxter equa- 
tions by exploiting the associativity of the extended Zamolodchikov-Faddeev 
algebrJl™OT 

S(9 12 )[I ® i^(0i)]S(0 12 )[I ® i^(fe)] = [I ® i^(# 2 )]S(<? 12 )[I ® i?^(0i)]5(0 12 ), (5) 

s(9 12 )[i®rP(9 1 )]s(9 12 )[i®t%(9 2 )} = R}{e x )®T^{e 2 ), (6) 

Sifix^T^)®^)] = [Ti(9 1 )^T^(9 2 )]S(9 12 ), (7) 

and 

Ri{9 l )®R}{6 2 ) = R}{e 1 )®Ri{6 2 ), (8) 

[Tf (0 2 ) ® I]5(^ 2 )[^(^) (g) I]S(0 12 ) = Tj(# 2 ) <g) Agfa), (9) 

[i®fP{9 2 )]s(e l2 )[i®Rl{e 1 )]s{e 12 )^Ri{e 1 )®fj{e 2 ), (10) 

[Tf(^)®I]5(^ 2 )[fj(^ 2 )®I] = [I®Tf(0 2 )]S(M[l®T^!)]. (11) 

We used here the convention (A <g) S)y = A^Bj for the tensor product and abbre- 
viated the rapidity sum 9\ 2 = 9\ + 9 2 and difference 9\ 2 = 9\ — 9 2 . Once again the 
same equations also hold for R <-> R and T «-> T. 

Apart from some discrepancies in the indices the equations 101-10 correspond 
to a more simplified, in the sense that there were no degrees of freedom in the 
d efec t and parity invariance is assumed, set of equations considered previously 
irP-i- 1 . For diagonal scattering it was argued m^ 1 that one can only have reflection 
and transmission simultaneously when S = ±1. IaP^la more general set up which 
includes all degrees of freedom was studied. A second set of equations l(%Jl-lfTT]l. 
which is not equivalent to ©-0 wa s found. It was shown that in the absence of 
degrees of freedom in the defect no theory which has a non-diagonal bulk scattering 
matrix admits simultaneous reflection and transmission. This result even holds 
for the completely general case including degrees of freedom in the defect upon 
a mild assumption on the commutativity of R and T in t hes e variables. It was 
further shown that besides S = ±1 also the Federbush modeP^and the generalized 
coupled Federbush model allow for R ^ and T ^ 0. However, when treating 
non-relativistic theories, the amplitutes depend only on the individual rapidities, 
such that one can fin d no n-trivial solutions with R ^ and T ^ even when the 
bulk theory is not free^U 

2.1.1 Multiple impurity systems 

The most interest situation in impurity systems arises when instead of a single 
one considers multiple defects, since that leads to the occurrence of resonance phe- 
nomena and when the number of defects tends to infinity even to band structures. 
Assuming that the distance between the defects is small in comparison to the length 
of the wire one can easily construct the transmission and reflection amplitudes of 
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the multiple defect system from the knowledge of the corresponding quantities in 
the single defect system. For instance for two defects one obtains 



T?(e)Tf(8) R a (e) = Ra (0) R?(9)T?(0)fm 



nap. 



fmf[{ff} Ra0{Q) = R0m RT(6)Tf(9)ff(9) 



These expressions allow for a direct intuitive understanding, for instance we note 
that the term [1 - R? (9)R? (0)}- 1 = T,n=i( R i ( )R?( d )) n simply results from the 
infinite number of reflections which we have in-between the two defects. This is of 
course well known from Fabry-Perot type devices of classical and quantum optics. 
For the case T = T, R = R the expressions Ijl2(l and coincide with the formulae 
proposed i When absorbing the space dependent phase factor into the defect 
matrices, the explicit example presented i for the free Fermion perturbed with 
the energy operator agree almost for T — T,R = R with the general formulae (1121) . 
They disagree in the sense that the equality of Rf^(0) and R^(9) does not hold 
for generic a, (3 as stated irl^l 

It is now straightforward to generalize the expressions for an arbitrary number 
of defects, say n, in a recursive manner 

rpotx ...a k ,Q\rpa k + i...U n ,n\ 

if (6) = K { ' 1 - jP—, 1< k < n, (14) 

Rf{0)=Rf 1 - ah (9) + ^ rP^- ^, Kk<n. (15) 

We encoded here the defect degrees of freedom into the vector a—{ai, ■ ■ ■ ,a n }. 
Similar expressions also hold for ff{0) = f * 1 - a "(9) and Rf{0) = £< n ~ a "(0). 

Alternat ivel y, we can define, in analogy to standard quantum mechanical meth- 
ods (see e.gP^J), a transmission matrix which takes the particle i from one side of 
the defect of type a to the other 

M i (e)-( ™ _1 -Rmmey^ m 

Jla{) -\-R«{-0)T?{-0)- 1 T?{-0)~ 1 )■ [W) 

Then alternatively to the recursive way l|14|l and l|15|l . we can also compute the 
multi-defect transmission and reflection amplitudes as 

Tf(0)= (f[Mi k (6)\ , Rf(6) = -(f[Mi k (9)\ (f[Mi k (6)\ . 

\k=l ) n \k=l / 12 \k=l J n 

(17) 

This formulation has the virtue that it is more suitable for numerical computations, 
since it just involves matrix multiplications rather than recurrence operations. In 
addition, it allows for an elegant analytical computation of the band structures for 
71 — ► oo, which we will however not comment upon any further. 
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2.2 Constraints from potential scattering theory 



As we argued above, in order to obtain a non-trivial conductance we are lead by 
integrability to consider free theories, possibly with some exotic statistics. Trying 
to be as close as possible to some realistic situation, i.e. electrons, we consider first 
the free Fermion, which, wi th a li ne of defect, was first treated irl^H The reafter it 
has also been considered irP^EZI arL( j3K| f rom different points of view. iJIIESOZI 
the defect line was taken to be of the form of the energy operator and i also a 
perturbation in form of a single Fermion has been considered. IrPZI we treated a 
much wider class of possible defects. 

Let us consider the Lagrangian density for a complex free Fermion -0 with I 
defects' 1 



(18) 



The defect is described here by the functions T> a ™ (-0, ip, dtxjj, dtip), which we assume 
to be linear in the Fermi fields ip,ip and their time derivatives. We can now pro- 
cee d in ana logy to standard quantum mechanical potential scattering theory (see 
and construct the amplitudes by adequate matching conditions on the 
field. We consider first a single defect at the origin which suffices, since multiple 
defect amplitudes can be constructed from the single defect ones, according to the 
arguments of the previous section. We decompose the fields of the bulk theory 
as %j)(x) = 0(x) ip+(x) + 0(— x) tj)-{x), with 0(x) being the Heavyside unit step 
function, and substitute this ansatz into the equations of motion. As a match- 
ing condition we read off the factors of the delta function and hence obtain the 
constraints 



il 1 (^+(x) - ifj-(x))\ x=0 



dV 



dip(x) 



d_ 

dt 



dV 



d(dti>(x)) 



(19) 



We then use for the left (— ) and right (+) parts of ip the well-known Fourier 
decomposition of the free field 



with the Weyl spinors 



(10 



(20) 



Uj {6) = -i-ySW) 



a We use the conventions: 

= (x ^ 1 ), = (in cosh 9, msinh 9), 



1 

1 



7 



1 

-1 



e -0/2 

e e/2 



7 



(2) 



(21) 



We adopt rclativistic units 1 = c = S = m K e 2 137 as mostly used in the particle physics 
context rather than atomic units l = e= /i=mRJ e/137 more natural in atomic physics. 
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and substitute them into the constraint 1)19(1 . Treating the equations obtained in this 
manner componentwise, stripping off the integrals, one can bring them thereafter 
into the form 

aj,-(0) = R^a^i-O) + T j (9)a j , + (9) , (22) 

which defines the reflection and transmission amplitudes in an obvious manner. 
When parity invariance is broken, the corresponding amplitudes from the right to 
the left do not have to be identical and we also have 

aj ,+(-9) = f.(9)a^(-9) + R0) aj , + (0) . (23) 

The creation and annihilation operators a\ (9) and (9) satisfy the usual fermionic 
anti-commutation relations {a,i(9\), aj(9 2 )} = 0, {a,i(9i), Oj(6 2 )} = ^SijS(9i 2 ). In 
this way one may construct the i?'s and T's for any concrete defect which is of 
the generic form as described in l|18|l . After the construction one may convince 
oneself that the expressions found this way indeed satisfy the consistency equations 
like unitarity QJ, O' an d crossing ©, (@J. Unfortunately the equations 
can not be employed for the construction, since they are not restrictive enough by 
themselves to determine the i?'s and T's. We consider now some concrete example: 

2.2.1 Impurities of Luttinger liquid type T>(ip, ip) = 4>(gi + g 2 ^f°)ip 
Luttinger 

liquidJSi 

are of great interest in condensed matter physics, which is one 
of the motivations for our concrete choice of the defect T>(ip, ip) = ip(gi + g 2 ~f°)ip. 
When taking the conformal limit of the defect one obtains an impurity which played 
a role in this context, see e.g!^, after eliminating the bosonic number counting 
operator. In the way outlined above, we compute the related transmission and 
reflection amplitudes 



- „ 4 Z (<7 2 + 5 icosh6»)e 2 ^ msinhe 
R 3 (9,9i,92,-y) = RM9x,92,y) = (4 + g2 _ 5 2 )sinh0 _ 4z(ffi+ff2COsh ^(24) 

R,(9,g u92 ,-y) = RM9U9„V) ^ (4 + g? „ g 2) sinM _ 4%1 _ g2 coshr ( 25 ) 



T j (6, gi , 92 ) = T j (9,g 1 , 92 )= {A + g 2_ g 2 )shlh9 _ il{gi+g2COshe) ■ (26) 



T 3 {6, gi ,g 2 ) = fj(9, gi ,g 2 ) = - ■ 2 [JZln L. TTT^TIT ■ ( 2 ?) 



(4 + 9% - 9 f) sinh6> 
(4 + g\ - gl) sinh 9 - Ai{g x - g 2 cosh 9) 

In the limit lim g2 ^o ^ (ip, ip) — ffiV'V': we recover the related results for the T/T's 
and R/R's for the energy defect operator. For this type of defect we present \T\ 2 
and |i?| 2 in figure 2 with varying parameters in order to illustrate some of the 
characteristics of these functions. 

Part (a) of figure 2 confirms the unitarity relation Q . Part (b) and (c) show the 
typical resonances of a double defect, which become stretched out and pronounced 
with respect to the energy when the distance becomes smaller and the coupling 
constant increases, respectively. Part (d) exhibits a general feature, that is when 
the number of defects is increased, for fixed distance between the outermost defects, 
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the resonances become more and more dense in that region such that one may speak 
of energy bands. 




Figure 2. (a) Single defect with varying coupling constant. |T| 2 and |i?| 2 correspond to curves 
starting at and 1 of the same line type, respectively, (b) Double defect with varying distance y . 
(c) Double defect with varying effective coupling constant B = arcsin(— 4gi/(4 + g\)). (d) Double 
defect = dotted line, eight defects = solid line. 



2.3 Laser fields interacting with quantum wires 

Let us now consider a more complex situation in which a three dimensional laser 
field hits the quantum wire polarized in such a wa y that it has a vector field com- 
ponent along the wire. Since the work of WeyP^, one knows that matter may be 
coupled to light by means of a local gauge transformation, which reflects itself in 
the usual minimal coupling prescription, i.e. dn — > d„ — ieA^, with being the 
vector gauge potential. The free Fermions in the wire are then described by the 
Lagrangian density 

C A = ${vfd» - m + erfA^ . (28) 
When the laser field is switched on, we can solve the equation of motion associated 

to (GEJ 

(iVfy - m + e^A^ = (29) 
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by a Gordon- Volkov type solution^! 

il>j(x,t). (30) 

Using now a linearly polarized laser field along the direction of the wire, the vector 
potential can typically be taken in the dipole approximation to be a superposition 
of monochromatic light with frequency u>, i.e. 

A(t) := Ax(t) = - f dsA (s) = -- [ dsE(s) = -— f cfe/(s) cos(cjs) (31) 
x Jo 2 J 2 J 

with f(t) being an arbitrary enveloping function equal to zero for t < and t > r, 
such that r denotes the pulse length. In the following we will always take f(t) = 
0(i)0(r — t) , with 0(x) being again the Heavyside unit step function. The second 
equality in H31(l . Aq(x, t) — xA(t), follows from the fact that we have to solve (|30|) . 

We want to comment on the validity of the dipole approximation in this context. 
It consists usually in neglecting the spatial dependence of the laser field, which 
is justified when xoj < c = 1, where a; is a representative scale of the problem 
considered. In the context of atomic physics this is typically the Bohr radius. In the 
problem investigated here, this approximation has to hold over the full spatial range 
in which the Fermion follows the electric field. We can estimate this classically, in 
which case the maximal amplitude is eE^/oj 2 and therefore the following constraint 
has to hold 

(^) 2 =4C/ P <1, (32) 

for the dipole approximation to be valid. Due to the fact that a; is a function 
of u>, we have now a lower bound on the frequency rather than an upper one as 
is more common in the context of atomic physics. We have also introduced here 
the ponderomotive energy U p for monochromatic light, that is the average kinetic 
energy transferred from the laser field to the electron in the wire. 

The solutions to the equations of motion of the free system and the one which 
includes the laser field are then related by a factor similar to the gauge transfor- 
mation from the length to the velocity gauge 

ipf(x, t) = exp [ixeA{t)\ ^(x) . (33) 

In an analogous fashion one may use the same minimal coupling procedure also to 
couple in addition the laser field to the defect. One has to invoke the equation of 
motion in order to carry this out. For convenience we assume now that the defect 
is linear in the fields ip and tp. The Lagrangian density for a complex free Fermion 
tp with I defects P a (^, ip, A^) of type a at the position x n subjected to a laser field 
then reads 

t 

C A D=C A + Y,^ an ^,4',^)S(x-x n ). (34) 

n=l 

Considering for simplicity first the case of a single defect situated at x — 0, the 
solution to the equation of motion resulting from (|34|l is taken to be of the form 
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ipf(x,t) = exp 



ie / dsAi(s,t) 



ipUx,t) = 



exp 



ie / dsAo(x, s) 



ipf~(x, t) = Q(x)ip^ + (x, t) + Q(—x)ip^_ (x, t) , which means as before we distinguish 
here by notation the solutions (|33[1 on the left and right of the defect, ipf_(x, t) and 
ipf+( x ,t), respectively. Proceeding as before, the matching condition reads now 

i A i i\ i A I +\\\ _ 9VAD(i',i>,A li ) 



i^Wft+foi) - ^_(x,t))\ x=0 - — 



(35) 

x=0 

It is clear, that in this case the transmission and reflection amplitudes will in 
addition to 9 and g also depend on the characteristic parameters of the laser field 

T(9, g, Eq, w, t) and R(6, g, E ,u,t) . (36) 

It is clear that the laser field can be used to control the conductance. For instance 
defects which have transmission amplitudes of the form as the solid line in figure 1 
(c), can be used as optically controllable switching devices. We can now turn to the 
central question which we address here, namely whether it is possible to generate 
harmonics in quantum wires. 



2-4 Analysis of the transmission amplitudes 

In order to answer that question, we first have to study the spectrum of frequencies 
which is filtered out by the defect while the laser pulse is non-zero. The Fourier 
transforms of the reflection and transmission probabilities provide exactly this in- 
formation 

1 r 

T(n,6,E ,u,T) = - / dt\T(9,E ,uj,t)\ 2 cos(nt), (37) 

T JO 

1 r 

K(n,0,E o ,w,T) = - / dt\R(d,E ,uj,t)\ 2 cos(nt). (38) 

T Jo 

When parity is preserved for the reflection and transmission amplitudes, that is for 
real defects with T>* — T>, we have |T| 2 + |i?| 2 = 1, and it suffices to consider T in 
the following. 



2.4-1 Type I defects 

Many features can be understood analytically. Taking the laser field in form of 
monochromatic light in the dipole approximation H31J) . we may naturally assume 
that the transmission probability for some particular defects can be expanded as 

00 

\T!{e,U p ,u:,t)\ 2 = ^ 2fe (0)(4(7 p ) fc sin 2fc M)- (39) 

k=0 

We shall refer to defects which admit such an expansion as "type I defects" . Assum- 
ing that the coefficients t2k{0) become at most 1, we have to restrict our attention 
to the regime 4/7 p < 1 in order for this expansion to be meaningful for all t. Note 
that this is no further limitation, since it is precisely the same constraint as already 
encountered for the validity of the dipole approximation i':>2l) . The functional de- 
pendence of H39fl will turn out to hold for various explicit defects considered below. 
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Based on this equation, we compute for such type of defect 

It is clear from this expression that type I defects will preferably let even multiples 
of the basic frequency lu pass, whose amplitudes will depend on the coefficients 
tik{&)' When we choose the pulse length to be integer cycles, i.e. r = 2irn/ui for 
n 6 Z, the expression in l|4()(l reduces even further. The values at even multiples of 
the basic frequency are simply 

Tj(2^e,u p ) = (-ir5>*(0) (u P ) k ( k _ ) , (41) 

fc=0 ^ ' 

which becomes independent of the pulse length r. Notice also that the dependence 
on Eq and to occurs in the combination of the ponderomotive energy U p . Further 
statements require the precise form of the coefficients t2k(9) and can only be made 
with regard to a more concrete form of the defect. 

2.4.2 Type II defects 

Clearly, not all defects are of the form (j39|l and we have to consider also expansions 
of the type 

\T H {e,E Q /e^,t)\ 2 ^ ]T ^ fc (0)-^-cos"M)sin 2fc (a;i). (42) 

k,p=0 

We shall refer to defects which admit such an expansion as "type II defects". In 
this case we obtain 

00 p 



_^2k+2p 



k,p=0 1=0 

( \ 

{2k + 2l)\tf k {9) {2k + 2iy4 P k +1 {6)E 



I J (-l) l + 1 TLU 2 + 2k 

\ 

(43) 



. ff [(2 9 ) 2 - (S) 2 ] "Wm - 1) 2 - (S) 2 ] , 

\q=0 " q=l I 

We observe from this expression that type II defects will filter out all multiples of 
u>. For the pulse being once again of integer cycle length, this reduces to 



^,^^)= E E M ^"? W ^ k+PE o p ( p i) ( 

h n — n n \ / \ 



k,p=0 1=0 

and 



2k + 21 

k + I — n 



(44) 



T n {(2n-l)uAEole)= £ (U,)** 

k,p=0 1=0 

'p\ (2/c + 2Z)!(2n- l)£ 2p+1 



I J {l + k-n + l)\(l + n + ky: 



(45) 
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which arc again independent of t. We observe that in this case we can not combine 
the Eq and uj into a U p . 



2.5 One particle approximation 



In spite of the fact that we are dealing with a quantum field theory, it is known that 
a one particle approximation to the Dirac equation is very useful and physically 
sensible when the e xte rnal forces vary only slowly on a scale of a few Compton 
wavelengths, see e.g@3 We may therefore define the spinor wavefunctions 



aUd) 



e -ipj-x 



2tT 2 P ° 



Uj (0) 



27T 2 p° 



(46) 



(47) 



With the help of these functions we obtain then for the defect system 

A AO)] 



,(x,t) :=Vi (*,*) 



Q(x)T*(6) hi, u ,g(x,t) + ^i, a ,- e (x,t)R*( 



®{-x) [^i, u ,e(x,t) + 9 ilU ,-o(x, t)R*(0)] 



(48) 



and the same function with u — + v. Since this expression resembles a free wave, 
it can not be normalized properly and we have to localize the wave in form of 
a wave packet by multiplying with an additional function, g(p,t) in H2()|) and its 
counterpart g(x, t) in Ij48(l . typically a Gaufiian. Then for the function &f u s (x, t) — 
g{x, t)^f u g (x, t), we can achieve that ||<fr|| = 1. 

2.6 Harmonic spectra 

We are now in the position to determine the emission spectrum for which we need 
to compute the absolute value of the Fourier transform of the dipole moment 



I dt ($f tUje (x,tfx&f tUte {x,t))expint 
Jo 



(49) 



We localize now the wave packet in a region much smaller than the classical estimate 
for the maximal amplitude the electron will acquire when following the laser field. 
We achieve this with a Gaufiian g(x,t) = exp(— x 2 /A), where A <C eE n /uj 2 . 

2. 7 An example: Impurity of energy operator type 

As mentioned this type of defect, i.e. T>(?p,ip) = g?pip(x) can be obtained in a limit 
from the defect discussed in section 2.2.1. Coupling the vector potential minimally 
to it yields 

V AD {j, V, = g${l + e/mrfA^ , (50) 
by invoking the equation of motion. 
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2. 7. 1 Transmission amplitudes 

We can now determine the reflection and transmission amplitudes as outlined above 



R % {0,g,A/e,y) = Ri(9,g,-A/e,-y) = g, A/e,-y) = R,(0,g,-A/e,y) = 

[yi-cosh 9]e- 2ly sinhe 
[1 - yicoshf?] - i\ [J + 1 + A 2 - y 2 A 2 } sinh(9 ' 



(51) 



We denoted the differentiation with respect to time by a dot. The transmission 
amplitudes turn out to be 

Ti(0, g, A/e, y) = f t (9, g, -A/e, -y) = T f (0, g, -A/e, y) = %{9, g, A/e, -y) = 



l~y 2 A 2 + {A- f) 2 



[1 - yAcosh9] ~ i[\ + 1 + A 2 - y 2 A 2 ] sinh6> 



sinh 9 

(52) 



Locating the defect at y — 0, the derivative of A does not appear anymore explicitly 
in (|51|) and (I52II . such that it is clear that this defect is of type I and admits an 
expansion of the form 1|39[) . With the explicit expressions (|51|l and l|52l) at hand, 
we can determine all the coefficients t2k(9) in H39fl analytically. For this purpose 
let us first bring the transmission amplitude into the more symmetric form 

, Tf o A i m2 M0,g)+a2(9,g)A 2 +a 4 (9,g)A i 

\Ti(9,g,Ae)\ = (53) 

ao(6»,g) +a 2 {9,g)A z +a i {9,g)A i 

with 

a ((9,g) = I6.9 2 + (4 + 5 2 ) 2 sinh 2 9, a (9, g) = {g 2 - A) 2 sinh 2 9, (54) 
a 2 (9,g) = 2g 2 (A + g 2 )smh 2 9, 0,4,(9, g) = g 4 sinh 2 6». (55) 

We can now expand |T(0, g, A)\ 2 in powers of the field A(t) and identify the co- 
efficients t2k(9,g) in l|39|) thereafter. To achieve this we simply have to carry out 
the series expansion of the denominator in (|53|l . The latter admits the following 
compact form 



J2c 2k (9,g)A 2k 1 (56) 



a {9,g)+a 2 (9,g)A 2 +a i {9,g)A 4 , n 

with c (9,g) = l/a (9,g) and 

, a \ C2k-2{9,g)a 2 (9,g) + c 2k -4(9,g)a 4 (9,g) 
C2k{9,g) = t^-^ , (57) 

for k > 0. We understand here that all coefficients c 2k with k < are vanishing, 
such that from this formula all the coefficients c 2k may be computed recursively. 
Hence, by comparing with the series expansion (|39|l . we find the following closed 
formula for the coefficients t 2k (9,g) 

t 2 k(0,g) = [a o (9,g)-a o (9,g)]c 2k (9,g) k > 0. (58) 
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The first coefficients then simply read 

to(e,g) = ^4 = \T(e, E() = o)\\ 



h(e,g) 



ao(6,g) 
(12(8, g) 



i-t {e, g )] = 



85 4 (4 + g 2 ) sinh 2 29 



a o (0,g) L ~ " u ^' a 'J (16g 2 + (4 + 5 2 ) 2 sinh 2 (9) 2 ' 
0-4(6, g) 0-2(8, g) 



02(8, g) a (8,g) 



t 2 (0,g), 



(59) 
(60) 
(61) 



and so on. It is now clear how to obtain also the higher terms analytically, but 
since they are rather cumbersome we do not report them here. 

2.7.2 Transmission amplitudes 

Having computed the coefficients t2k, we can evaluate the series (|40f) and (|41H in 
principle to any desired order. For some concrete values of the laser and defect 
parameters the results of our evaluations are depicted in figure 3. 




15 PVC020 



2S 30 35 40 



Figure 3. Fourier transform of the transmission probability for a single (a) and double (b) defect 
with Eq = 2.0, g = 3.5, 8 = 1.2, w = 0.2. Harmonic emission spectrum for a single (c) and double 
(d) defect with E = 2.0, g = 3.5, 9 = 1.2, uj = 0.2, A = 6. 



The main observation from part (a) is that the defect acts as a filter selecting higher 
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harmonics of even order of the laser frequency. Furthermore, from the zoom of the 
peak regions, we see that there are satellite peaks appearing near the main harmon- 
ics. They reduce their intensity when r is increased, such that with longer pulse 
length the harmonics become more and more pronounced. Wo also investigated 
that for different frequencies u> the general structure will not change. Increasing 
the field amplitude E , simply lifts up the whole plot without altering very much 
its overall structure. We support these findings in two alternative ways, either by 
computing directly l|37|) numerically or, more instructively, by evaluating the sums 
|@DJ and g]}. 

Part (b) shows the analysis for a double defect system with one defect situated 
at x = and the other at x = y. The double defect amplitudes are computed 
directly from (|12Jl and (|13[1 with the expression for the single defect (|51|> and (|52(l . 
Since now both A and A appear explicitly in the formulae for the i?'s and T's, it 
is clear that the expansion of the double defect can not be of type I, but it turns 
out to be of type II, i.e. of the form l|42(l . Hence, we will now expect that besides 
the even also the odd multiples of u> will be filtered out, which is indeed visible in 
part (b) for various distances. Here we have only plotted a continuous spectrum for 
y = 0.5, whereas for reasons of clarity, we only drew the enveloping function which 
connects the maxima of the harmonics for the remaining distances. We observe 
that now not only odd multiples of the frequency emerge in addition to the ones in 
(a) as harmonics, but also that we obtain much higher harmonics and the cut-off is 
shifted further to the ultraviolet. Furthermore, we observe a regular pattern in the 
enveloping function, which appears to be independent of y. Similar patterns were 
observed before in the literature, as for instance in t he c ontext of atomic physics 
described by a Klein-Gordon formalism (see fig ure 2 irPSj). 

Coming now to the main point of our analysis we would like to see how this 
structure is reflected in the harmonic spectra. The result of the evaluation of 
is depicted in figure 3 parts (c) and (d). We observe a very similar spectrum as 
we have already computed for the Fourier transform of the transmission amplitude, 
which is not entirely surprising with regard to the expression (|49|l . The cut-off 
frequencies are essentially identical. From the comparison between X and the 
enveloping function for T we deduce, that the term involving the transmission 
amplitude clearly dominates the spectrum. 

The important general deduction from these computations is of course that 
harmonics of higher order do emerge in the emission spectrum of impurity systems, 
such that harmonics can be generated from solid state devices. 



2.7.3 Relativistic versus non-relativistic regime 

In the previous sections we have been working in an intensity regime which is clos e to 
the damage threshold of a solid, according to the experimental observations i 
This allowed us to see the maximum effect with regard to harmonic generation 
which at present might be visible from experiments. However, it is also interesting 
to investigate situations which are not experimentally feasible at present and of 
course lower intensity regimes. 

In order to judge in which regime we are working and whether there are rel- 



Proceeding 6-th international workshop on CFT and IM (Chernogolovka, 2002) 15 



ativistic effects, let us carry out various limits. First of all we recall a standard 
estimation according to which the relativistic kinetic energy is close to the classical 
one when one is dealing with velocities v 2 <C 3/4c 2 . This is the same as saying 
that the kinetic energy is much smaller than the rest mass -Ekin "C tuqc 2 . Making 
now a rough estimation for the system under consideration, we assume that the 
total kinetic energy is the one obtained from the laser field, i.e. the ponderomotive 
energy U p . We also ignore for this estimation any sophisticated corrections, such 
as possible Doppler shifts in the frequency, etc. Then the non-relativistic regime is 
characterized by the condition U p <C 1. 

Based on our previous observation that T and X exhibit a very similar be- 
haviour, it will be sufficient here to study only the T in the different regimes, 
which will be easier than an investigation of the full emission spectrum l(49jl • From 
our analytic expression l|4*T|) . we see that for a type I defect the quantity 7} becomes 
a function of U p , such that the regime will be the same when we rescale simultane- 
ously Eq and lu. Accordingly we evaluate numerically the Fourier transform Qllfl. or 
equivalently compare against our analytical expression l|41|) , and depict our results 
in figure 4a. 




Figure 4. Absolute value squared of the Fourier transform of the transmission probability, (a) 
Single defect for various values of U p with g = 3.5, 9 = 1.2. (b) Double defect for varying values 
of E with g = 3.5, d = 1.2, v = 0.2. 

We observe that when passing more and more towards the relativistic regime the 
cut-off is increased. The other feature one recognizes is that the modulating struc- 
ture in the enveloping function of the harmonics becomes more pronounced. One 
should also note, in regard to l|32|) . that the multipole structures might become 
more and more important in the relativistic regime. 

Let us now perform a similar computation for the double defect. From the 
expressions (144(1 and l(45|l we see that now 7}/ is not just a simple function of U p 
and therefore even being in the same regime the behaviour will be different when Eq 
and io are rescaled. We alter in that case the regimes by rescaling Eq and keeping 
the frequency fixed. Our results are depicted in figure 4b. 

Similar as for the single defect we see that the cut-off is increased and the 
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modulating structure in the enveloping function becomes more emphasized when 
we move towards the relativistic regime. In addition we note that the difference 
between the even and odd harmonic becomes larger with increasing U p . This effect 
is more extreme for the low order harmonics. 

As a general observation we state that there are not any effects which seem to 
be special to the relativistic regime, but the transition to that regime seems to be 
rather smooth. 

We will now consider a quite different application, which allows, however, to use 
integrability not only in a constraining but in a completely constructive manner. 



3 Fractional quantum Hall systems 

The quantum^! and in particular the fractional! quantum Hall effect have at- 
tracted an enormous amount of attention b oth from theorisll^SI and experimental- 
ists (for some very recent experiments see e.gl^j). The key observation is that when 
subjecting an electron gas confined to two space dimensions to a strong uniform 
magnetic field, the transverse (Hall) conductance takes on preferably certain char- 
acteristic values G — e 2 /hv, whereas the longitudinal conductance vanishes at these 
plateaux in complete analogy to the classical Hall effecl)^! The filling fractions v 
are distinct universal, in the sense that they are independent of the geometry or 
type of the material, rational numbers, which can be determined experimentally 
to an extremely high precision. Many, but not all, of the experimentally observed 
filling fractions are part of Jain's famous sequence 

(se<P 

and references therein) 

TTL 

v= — m,p/2 = 1,2,3,... (62) 

mp ± 1 

which results as a theoretical prediction from a composite Fermion theory. 

In the following we will show that these universal numbers also quantize the con- 
ductance of quantum wires when described by minimal affine Toda field theorie^^l. 



3.1 Conductance in the high temperature regime 

There exist two established theoretical descriptions to compute the conductance. 
The first is based on a linear response theory in which one essentially needs the 
Fourier transform o f th e current-current two-point correlation function. _ This so- 
called Kubo formulcPEI has been adopted to a situation with a boundary^! How- 
ever, since this set-up only captures effects coming from the constriction of the wire 
a generalization which includes defects was needed, which we proposed i as 



G a (l/T) = - lim — - / dt e*«" (J(t)Z a J(0)) T ,m ■ (63) 

u->0 2ldlT A J_ 00 

Here the defect operator Z a enters in-between the two local currents J within 
the temperature T and mass to dependent correlation function. The Matsubara 
frequency is denoted by ui. 

The other possibility of determining the conductance is a generalization of the 
Landauer-Biittinger transport theory picture. Within this framework a proposal 
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for the condu ctance through a quantum wire with a defect (impurity) has been 
made irPHEOI 

G Q (1/T) = ]T ton | f°° de\ P W,T,^\T?(e)\ 2 -pU9,T,ti)\f«(e) 
i (p'-mD^o 2 J_ QO L 

which we only modify to accommodate parity breaking 5 ^. This means we allow the 
transmission amplitudes for a particle of type i with charge qi passing with rapidity 
9 through a defect of type a from the left T" (9) and right Tf (9) to be different. 
The density distribution function p r i {9,T 1 /Zj) depends on the temperature T, and 
the potential at the left p, 1 - and right p\ constriction of the wire. 

In both descriptions (|63|l and (|64|l one can employ non-perturbative methods 
of integrable models, leading in (|63|l to an exac t expr ession for the current-current 
correlation functions (. . )T,m from a form factoJ ' ' • expansion and in (|64l) for the 
density distributions pt from a thermodynamic Bethe ansatJ*^ analysis. 

Here we will not consider ()63(l and only present a more detailed study of H64JI 
when impurities are absent and parity is preserved. In order to compute G in a one 
dimensional quantum wire, we simply have to determine the difference of the static 
charge distribution at the left and right constriction of the wire, which we assume to 
be at the potentials p\ and respectively. Then, to obtain the direct current Jj for 
each particle of type i with charge qi, we have to integrate the density distribution 
functions p r i {9,T,p i ) of occupied states over the full range of the rapidities 9 and 
the conductance simply reads 

G(l/T) = V ton — !— Z;(l/T, Apt = (65) 



= E A |™ o gf- / d0 T > /4) - Pi(0, T, p\)} . (66) 

Hence, the main task in this approach is to determine the density distribution 
functions p r i {9 1 T, p{) of occupied states. It is remarkable that in the context of 
integrable models, despite the fact that these functions are neither Fermi-Dirac 
nor Bose-Einstein, there exist approaches in which they can be computed non- 
perturbatively, i.e. the thermodynamic Bethe ansata^l 

We briefly recall how this works. The central equations of the TBA relate the 
total density of available states pi(9,r) for particles of type i with mass m, as 
a function of the inverse temperature r = l/T to the density of occupied states 
Pl(9,r) 

Pi (9, r) = ^ coshtf + * p]]{9) . (67) 

3 

By (/ * g) (9) := l/(2vr) / d9'f(9 - 9')g(9') we denote as usual the convolution of 
two functions. There are only two inputs into the entire TBA analysis: first the 
dynamical interaction, which enters via the logarithmic derivative of the scattering 
matrix tpij{9) = —id\nSij{9)/d9 and an assumption on the statistical interaction 
gij amongst the particles i and j on which we comment further below. For the 
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moment we chose this interaction to be of fermionic type. The mutual ratio of 
the two types of densities serves as the definition of the so-called pseudo-energies 
£i(0,r) 

Pi (6,r) 1 + [ ' 

which have to be positive and real. At thermodynamic equilibrium they can be 
computed from the non-linear integral equations 

rim cosh0 = £i(0,r, fa) + rfa +^[^3 * + e~ e *)](6) , (69) 



where r = m/T, mi — > m;/m, pi — ► Pi/m, with m being the mass of the lightest 
particle in the model and chemical potential fa < 1. As pointed out already i 
(here just with the small modification of a chemical potential), the comparison 
between ((69(1 and 167(1 leads to the useful relation 



Pi(6,r,fa) = 



1 / cfei(6>, r, /ii) 



2tt 



(TO) 



Here e(0) = 0(6>) - 8(-0) is the unit step function, i.e. e(0) = 1 for > and 
e(0) = — 1 for < 0. In equation ((68(1 . we assume that in the large rapidity regime 
the density p^(0,r,fa) is dominated by the last expression in 170(1 and in the small 
rapidity regime by the Fermi distribution function. Therefore, from ((68(1 follows 



p^(9,r,fa) = 



-£i(6,r,fii) 



1 -)- e -£i(8,r,fii) 



pi(9,r, p^ 



^ e ( )4 ln [1 + exp(-£;(0, r, pi))] 
2irr do 



(71) 
(72) 



Using this expression in equation 1(66(1 . we can approximate the direct current in 
the ultraviolet by 



YanUr, A fa) ~ - — 
1 — >o 47rr 



dOhi 



1 + exp(-e j (6»,r, p[)) 



l + exp(- £i (0,^/4)) 



de(fl) 
d(9 



(73) 



after a partial integration. Taking now the potentials at the end of the wire to be 
p\ = —pb\ = pi/2 we carry out the limit Api — > in ((661) with the help 1'HospitaPs 
rule and the conductance becomes 



l imGi (r)~£- 
r^o Zirr 



dO 



1 



d£i(e,r, pi/2) 



Noting that de(9)/d9 = 25{6), we obtain 

q t 1 



1 + exp[£i(0,r,O)] dpi 

efei(0,r, /ij/2) 



lim GAr 1 
r^a irr 1 + expej(0,r, 0) 



Mi=0 



de(fl) 
d6> 



The derivative cfei(0,r, pi/2)/dpt can be obtained by solving 



cfe t (0,r, /i,/2) 
dp k 



de J (0,r,p j /2) 



exp£j(0,r, pi/2)} 



dp k 



(74) 



(75) 



(76) 
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N l3 = — lim [ln[^- i-OySij (6)}} . (77) 



which results from performing a constant TBA analysis on the «fc-derivative of 16911 
in the spirit oP. At this point only the asymptotic phases of the scattering matrix 
enter via 

1 

2ni e^+c 

In principle we have now all quantities needed to compute the conductance, but 
to solv e (|76f) for the derivatives of the pseudo-energies is somewhat cumbersome, 
ggJSZEl f or such a computation. Nonetheless, we can elaborate more on equation 
(|76[) and simplify the procedure further. For this purpose we introduce the quantity 

such that we can re-write equation (|76() equivalently as 

MijY jk = with My := N l3 - (1 + e £i )<^ (79) 

where the pseudoenergies satisfy the constant TBA equations 

e- ei = JJ(i + e - e O Ww ■ (80) 

3 

Returning now to dimensional variables, i.e. replacing l/2w — » e 2 /h, the conduc- 
tance at high temperature in terms of the filling fraction v then simply results 
to 

e 2 1 , 

G(0) = —v with v = ^£ J <li{M- 1 ) ij . (81) 

id 

This means we have reduced the entire problem to compute filling fractions simply 
to the task of finding and inverting the matrix M. This is done in two steps: 
First from the asymptotic phases of the scattering matrix we compute Nij and 
subsequently we solve the constant TBA equations JSUJ). Then it is a simple matter 
of inverting the matrix (|79|l and performing the sums in (|81|l . 

In the context of the fractional quantum Hall effect one encounters very often 
particles which obey some exotic (anyonic) statistics. So far we have assumed our 
particles to obey fermionic type statistics as this choice is most natural for the 
investigated theorieJ^. However, one c an e asily implement more general statistics 
by adding a matrix gij to the iV-matriJ^I. 

The formula l|81(l reminds of course o n the w ell-known expressions for th e con- 
ductance as may be found for instance iJSSEH In that context it was found^EsO 
that Jain's sequence l|62l) can be obtained simply from the (m x m)-matrix 

M, . l> : <>, ., • (82) 

For this we have to take qi — 2 V i in our expression l|81|l . We will now demonstrate 
that the sequence (16211 can also be obtained in a more surprising way from fairly 
complicated matrices, even with non-rational entries, which result directly in the 
way indicated above, namely from a TBA analysis of minimal affine Toda field 
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theoriei^l Each Toda theory is associated to a Lie algebra g of rank £ and it is 
well knowrJSH that in that case N is an (I x ^)-matrix which is of the general form 

N tJ = S tJ - 2(K-% , (83) 

where K s is the Cartan matrix related to g ( see e.g. The solutions to the 

constant TBA equations are also knowrPSlaSl for most cases. In the ultraviolet 
limit these theories possess Virasoro central charge c = 21/ (H + 2), with H being 
the Coxeter number of the Lie algebra g. 

3.2 Fractional filling fractions from minimal affine Toda field theory 

Here we only reproduce the most important and stable subsequence of Jain's hier- 
archy, that is Ij62(l for p = 2. Let us start with some concrete examples to illustrate 
the working of our formulae. The first member of this series, that is filling fraction 
v = 1/3 is one of the best studied example a I I lof the fractional quantum Hall 
effect. We show now that this particular filling fraction results in the ultraviolet 
limit from an A3-affine Toda field theory. Specializing the general expression l|83|l 
to the ^3-case, the solutions to the constant TBA equations (18011 are simply 



(84) 



Then, the inverse of the M-matrix 



Mij = Sa - 2(#7% - <%(1 + e e * ) (85) 



is computed to 




(86) 



From the fact that the ^-mi nima l affine Toda field theories can also be viewed 
as complex sine-Gordon modelgSSl we knowED that the charges in this theory are 
q% = qs = 1, #2 = 2, such that l(HT|) yields 

va 3 = 1/3 . (87) 

This is not entirely surprising as it is knowrP^I that the fractional quantum Hall 
effect with this filling fraction can be described successfully in terms of a c = 1 
conformal field theory (CFT), as in the case at hand. The next example, i.e. As- 
minimal affine Toda field theory, yields a less expected answer, even more since 
the M-matrix contains non -ration al entries. With l|83|) for A5 the solutions to the 
constant TBA equations arJ^HEHI 

e- 1 = e £5 = 1 + V2, e £2 =e ei =2 + 2V2, e £3 = 3 + 2\[2 . (88) 

Assembling this into the M-matrix, it is clear that it will contain non-rational 
entries. Evidently this matrix is not of the form l|82(l and certainly falls out of 
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the classification scheme based on integral latticed^. Nonetheless, it will lead to a 
distinct rational value for v. We compute the inverse of M to 



(( 



M~ 



6^2) 
11 



( 




6 



17 

2 V2 



(^H ( 



2 

7-5 \/2 
4 



(6V2- 



15- 



i/5 



11 



(89) 



Remarkably when taking into account thaJ^gi = q§ = 1, q% = q^ = 2, q$ 
obtain by evaluating Ij81|l for the matrix (|89|l the simple ratio 

va, = 3/8 • 



= 3, we 



(90) 



We will now turn to the gen eric ca se. Taking the general solutions of the constant 
TBA equations into accountpSESl a nd using a generic expression for the inverse 
of the Cartan matrix K^j — min(i, j) — ij/(£ + 1) in Ij83(l . the M-matrix for an 
^2^+1 -minimal affine Toda field theory can be written generically as 



Ma = - 

3 £ 



ij 



- - 2min(i, j) - 5 ir 



(2^+4) 



(» + 2)7T 

2£+4 



sm 



2£+4 



(91) 



As already indicated by the previous example this matrix is not of the form 
and does not fit into the classification scheme proposed irP^l. According tcE J v. < - 
have the charges 



q2£+2- 



and 



qi = i for i < £ + 1 



(92) 



As can be guessed from (|89f) . it is not evident how to express the inverse in terms 
of a simple closed expression. We can, however, invert l|91|l case-by-case up to very 
high rank and we obtain from (|81() together with l|92|) the sequence 



VA 2 p 4. 1 — 



l+l 
2£ + A 



(93) 



Taking now £ — 2m — 1, we obtain as a subsequence of this the most stable part of 
Jain's sequence Ijfi2|l with p = 2 



VA — = 2^TT ' (94) 

In summary: The conductance of a quantum wire which is described by a massive 
A2£+i-minimal affine Toda field theory possesses in the high temperature regime, 
in which the model turns into a conformal field theory with Virasoro central c = 
(fit + 1) / (£ + 2), a filling fraction equal to ^9fy) . In particular for £ = 2m — 1, we 
obtain the Jain sequence \94\j - 

In order to illustrate the main idea we just presented here the most stable of 
the Jain sequence. To see how to obtain other sequences we refer the reader tcPl 
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4 Conclusions 



We have presented two concrete applications in which non-perturbative methods 
developed in the context of integrable quantum field theories can be used to evaluate 
physical quantities. We predict the generation of harmonic spectra from solid state 
devices and show that filling fractions occurring in fractional quantum Hall systems 
can be obtained from minimal afhnc Toda field theories. 
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